SUNY-NTG-92/40 



Mesonic Correlation Functions in the Random Instanton Vacuum 

E.V. Shuryak and J.J.M. Verbaarschot 

Department of Physics 
SUNY, Stony Brook, New York 11794 

Abstract 

A general model-independent discussion of mesonic correlation functions is given. We 
derive new inequalities, including one stronger than Weingarten's inequality. Mesonic 
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in all channels, but in several channels the model works with amazing accuracy up to 
distances of 1.5 fm. 
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1. Introduction 



1.1. Motivations and outline of the paper 

Correlation functions are a major tool for understanding the structure of any type 
of matter, solid or liquid. The QCD vacuum, being a complicated ensemble of interact- 
ing quark and gluon fields, has also been studied this way in the framework of various 
theoretical approaches based directly on quantum field theory, especially, by numerical 
simulations on a lattice. However, in the past, those studies were mainly focused on 
the long-distance behaviour of these functions, which is determined by the masses of the 
lowest hadronic states. Only recently was it realized, that extremely important dynam- 
ical information is contained in the behaviour of these functions at smaller distances of 
a; ~ 0.1 - 1.0 fm. 

Another explanation of why these point-to-point correlation functions are so funda- 
mental for QCD goes as follows. Due to confinement, one cannot directly study quark- 
quark scattering, but the same information is contained in the correlation functions. At 
the same time, these functions can also be obtained by inserting a set of physical (asymp- 
totic) states and using dispersion integrals. In this way, one can learn a lot about the 
space-time structure of the effective interquark interaction. The available phenomenolog- 
ical information on correlation functions, as well as the current theoretical situation, was 
recently reviewed by one of us [Q. 

Although some questions related to the short-distance behaviour of the correlators 
can be traced back to 'current algebra', the Wilson operator product expansion (OPE) 
and other developments in the sixties, a systematic study of the relevant non-perturbative 
QCD dynamics was initiated by the so called QCD sum rules 0. Using the OPE as a 
practical tool to extract information on correlators at small distances, these authors have 
found that it works nicely in some cases, but fails in others. 

In short, the general lesson one can draw from these studies is that the short range 
behaviour of the correlation functions is qualitatively different for different channels. 
Non-perturbative effects cause deviations from the perturbative behaviour, that differ 
drastically in sign and magnitude for different channels. The resulting picture obviously 
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contradicts to any simplistic model of strong interactions. In particular, universal color 
confinement is far from being the only non-perturbative force between quarks and anti- 
quarks. 

The main objective of the present work is to investigate the role of instanton-generated 
effects. By now, it is well known, that at least some of their qualitative features (briefly 
reviewed in section 1.6) are in good correspondence with empirical observations for the 
scalar and pseudo-scalar channels. Moreover, numerical studies in the framework of the 
'interacting instanton approximation' (IIA) show very promising results, outperform- 
ing the QCD sum rules even for the vector and axial correlators. In this work, we report on 
much more quantitative results for many of these functions and establish more definitely 
the limitations of this approach. 

Another motivation for these studies is that one can get a similar set of results from 
lattice simulations. The first work in this direction has been reported recently Q]. In 
the last section of this work, the comparison of the results obtained within these two 
approaches shows that they are surprisingly consistent with each other. 

Many details of our approach have already been explained in the first paper of this 
series so we do not need to repeat them here. Let us only remind, that we use the 
simplest possible model, the 'random instanton vacuum' (RIV), as a kind of 'benchmark' 
for future investigations of instanton dynamics. The main result of this first paper is the 
derivation of a practical approximate formula for the quark propagator. More specific re- 
sults for the averaged propagator can, roughly speaking, be interpreted as the appearance 
of some 'quark effective mass', which happens to be about 300 MeV in agreement with 
the old 'constituent quark model'. However, as we will show below, it does not mean 
that this model provides a reasonable approximation to correlation functions. The quark 
propagators contain 'hidden components' ! They are not seen in the average propagators, 
but they show up if the square (for mesonic correlators) or the cube (for baryonic corre- 
lators) of the propagators is averaged over the gauge field configurations. In other words, 
we have found a strong interaction between quarks, which is not small as compared with 
the 'constituent masses' and is strongly channel dependent. 

The paper is organized as follows. In section 1.2 we remind some general 'kinematical' 
properties of the correlators and introduce the currents and the terminology used. Then, 
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in section 1.3, we specify how the various correlators are connected to the propagator 
decomposition into spin-color components. These formulae can be used, e.g., to make the 
simplest Vacuum dominance approximation', which allows us to get at least the correct 
sign of the non-perturbative corrections to various correlators induced by chiral symmetry 
breaking. In this subsection we also derive some general inequalities between hadronic 
correlators, including one which is stronger than Weingarten's inequality 0. In section 
1.4, we give schematic spectral functions for the various flavored channels. A general 
discussion of unflavored correlators is given in section 1.5. We finish the introductory 
part of the paper by reminding the reader of the qualitative features of the instanton- 
induced effects on correlators (see subsection 1.6). 

The main body of the paper is essentially a one-by-one discussion of our results for 
the various channels, starting with a rather complete set of 'flavored' currents (section 2) 
supplemented by the most interesting 'unflavored' channels (section 3). In all cases our 
results are compared with phenomenological expectations (based on a detailed review []I|) 
or, if they are absent, with the available predictions from QCD sum rules. The comparison 
with lattice data and a general discussion of the results (including specific fits of masses 
and coupling constants for various mesons) can be found in section 4. 



1.2. General properties of the correlators 

First of all, correlators of two currents can be divided into diagonal and non-diagonal 
correlators, depending on whether they involve the same currents or two different currents. 
It is also convenient to recognize two different types of currents, which we call 'flavored' 
and 'unflavored' currents. The correlation function of the flavored currents {e.g. ud or us, 
but also uu — dd) involves only diagrams in which both quark and anti-quark propagate 
between x and y. The correlation function of the unflavored currents (like uu) also 
receive contributions from two-loop diagrams, in which one quark line goes from x to x 
and another quark line goes from y to y. The one-loop case, which is much simpler {e.g. 
two- loop correlators are considered rarely in lattice QCD calculations), will be studied 
first. In line with this, we differ from the usual SU (3) notation and organization in which 
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the correlators are classified according to their S'C/(3)— representation. Therefore, the 77 
and 77' mesons are discussed separately (see sections 1.5 and 3.1) rather than together 
with the 71 and K correlators. Another important channel in which 'two-loop diagrams' 
are crucial is the iso-scalar scalar a channel (see sections 1.5 and 3.2). 
In this paper we will consider the flavored currents 



= ud, (1.1a) 

= ui'f^d, (1.1b) 

Ji = M7^75C?, (1-ld) 

^ = ua^^^d, (l.le) 



and their counterpart with the u— quark interchanged by an s— quark. They obey the 
hermiticity relation 

(dVu)^ = uTd. (1.2) 

The study of unflavored currents will be restricted to the following (Hermitian) cur- 
rents: 



r 


— —i={uu + dd), 
v2 


(1.3a) 




= —= {ui^5U + di^sd + si^5s) , 
v3 


(1.3b) 


f 


= —7= {uij5U + di'j^d - 2si75s) , 
V6 


(1.3c) 


J 1^ 


= -^{u-fi,u + d'yi,d), 


(1.3d) 
(1.3e) 
(1.3f) 



The two-loop diagrams contributing to the last three currents will be ignored (Zweig's 
rule), and therefore they will be discussed together with the flavored currents. 
The correlators are defined by the time-ordered product 

Il{x)=<0\Tu{x)jimO>, (1.4) 



where the Dirac indices are suppressed. If we use completeness, and PCT and transla- 
tion invariance of the vacuum state, the correlator can be written in terms of a spectral 
representation as 

n(a;) = / ^(^(a;o)e-^- + e{-Xo)e^'i^)5{q' - a')p{a^)da\ (1.5) 
where p(cr^) is the spectral function defined by 

p{<j^) = (27r)3$:5V-g„)<0|j^(0)|n><n|j^(0)|0>, (1.6) 

n 

with Dirac indices as given by the two currents. The integral over d^q yields a Feynman 
propagator of mass — cr^, which for spacelike distances is given by 

D{a,x)^-^K^{ax), (1.7) 

where Ki[ax) is a modified Bessel function. The integral in (1.5) can be rewritten as 

/■oo 

n(x)=/ da'^D(a,x)p(a'^). (1.8) 

^0 

For spacehke x, the r.h.s of (1.8) can be analytically continued to Euclidean spacetime 
without hitting any singularities. This equation can thus be used to relate the spectral 
function defined in Minkowsky space-time to a correlator in Euclidean space-time. For 
later use we also give its Euclidean Fourier transform: 

[ d^xe'^^U(x) = r ds-$^. (1.9) 
J Jo + s 

On the other hand, again using PCT invariance, the spectral function can be written 
in terms of the discontinuity of the correlator in Minkowsky space-time 

27rp(q^) = 2Imi J < 0\TjA(x)jUo)\0 > ■ (1-10) 

The Euchdean Fourier transform of (1.8) with the spectral density replaced by the inte- 
gral in (1.10) yields the standard dispersion relation, in which the correlator evaluated 
at Euclidean momenta is given by an integral transform of the imaginary part of the 
correlator in Minkowsky space time. 

Finally, let us consider the matrix element of the current in (1.6) in more detail. For 
the currents defined in (1.1) obeying the Hermiticity relation (1.2) and the hermitian 
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currents (1.3) the relations 

< 0\dTu\n >*=< 0\u-f5T-f5d\n >, (1.11a) 

< Q\uTu\n >*=< 0|M75r75M|n > . (1.11b) 

follow from PCT invariance. Therefore, the matrix elements of a Hermitian current are 
real if the gamma matrix structure is spin flipping, i.e. T = S, P or T, and they are 
imaginary if the gamma matrix is spin nonflipping, i.e., T = V oi A. 

1.3. Propagator decomposition and inequalities 

In the previous paper we already mentioned, that the quark propagator contains 
both chirality-fiipping and non-fiipping components. Some of them survive the averaging 
over configurations, whereas others, called the 'hidden' components, are nonzero only if 
a square (or cube) of the propagator is averaged. The 'hidden' components generate the 
effective interaction between quarks. 

Thus, the study of mesonic or baryonic correlators actually reveals these 'hidden com- 
ponents'. Generally speaking, they are by no means small as compared with the 'visible' 
components. In other words, the interquark interaction is generally of the same order of 
magnitude as the 'constituent quark masses' studied in 0. 

In this section we specify the 'kinematics' of this statement for various channels, and 
briefly consider the relevant dynamics. The correlators, evaluated in Euclidean space-time, 
can be expressed in terms of the Euclidean propagator S^l{x, y). The general structure of 
the correlator is 11 =< TT{S{x,y)rS{y,x)r) >, where the average is over all gauge field 
configurations. Specifically, we discuss the following diagonal flavored correlators: 



Us = <TT{S{x,y)S{y,x))>, (1.12a) 

Up = <Tr(5(x,y)z755(y,x)275) >, (1.12b) 

Ul^ = <TTiS{x,y)^,S{y,xh,)>, (1.12c) 

n;^^ = <TT{S{x,y)-f,-f,S{y,x)j^-f,)>, (1.12d) 
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(1.12e) 



As we will see below, only two non-diagonal flavored correlators are non-zero. They are 
given by: 

nj"^ = <Tr(5(x,y)z755(l/,a:)7M75)>, (1.13a) 

Ur = <MS{x,y)j,S{y,x)a^,> . (1.13b) 



The propagator can be decomposed as 

5jt = E«fr,,., (1.14) 

i 

where the Fj are a complete set of 16 Hermitian gamma-matrices given by, 

Ti = 1, 75, 7^, ^7^75, and ^7^7^ (fi ^ v). (1.15) 

The coefficients a, are ^[/(S)— color matrices. Before rewriting the correlators in terms 
of the ttj we point out a relation given by Weingarten for the Euclidean propagator in 
backward direction 

S{,x,y) = --i^S\y,x)-ir^, (1.16) 

which allows us to express the diagonal correlators defined above in terms of the positive 
definite quantities 

kl 

In the V , A and T case the Dirac sum over all Qi belonging to a specific structure, is 
denoted by Qy, Qa and Qt, respectively. 

The pseudoscalar (vr) correlator, 11'^, is simply given by the sum of all coefficients 
squared 

= Qs + Qp + Qv + Qa + Qt- (1.18) 

In the scalar (S) case, the correlator II"^ is given by the difference of the spin non-fiipping 
and the spin flipping components: 

U'^ = -Qs-Qp + Qv + Qa-Qt- (1.19) 
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As a result, we find Weingarten's inequality that the pseudoscalar correlator should 
exceed the scalar one. From the experimental values of the pion and the scalar meson 
masses it is clear that this inequality should be satisfied at large distances. What is 
nontrivial is that it holds at all distances. One expects that for x > 0.5 fm the scalar 
correlator is practically zero, which implies there should be a very precise compensation 
between different components of the propagator. 

The vector (p) correlator, H^^, and the axial (ai) correlator, H^^, are given by 

-Ul^ = +AQs-AQp + 2Qv-2Qa, (1.20a) 
-n;^^ = ■~AQs + ^Qp + 2Qv-2Qa. (1.20b) 



At small distances both correlators are dominated by the contribution of the free propa- 
gator and are very similar. This is no longer the case at distances of the order of 1 fm |l|] 
which implies a fine tuning of the components of the propagator. The tensor correlator, 
given by 

K^^u = -QQs - QQp + 2Qt, (1.21) 

does not receive a contribution from the free propagator and is less singular at short 
distances than the four other correlators. The above relations can be inverted, which 
allows us to express the components Qi in terms of the mesonic correlation functions. Only 
the pion correlator always enters with a positive coefficient which implies that, if there is 
only one Goldstone boson, it necessary should be the pion. This is in agreement with the 
Vafa-Witten theorem [^] that vector symmetries cannot be spontaneously broken. 
Note that these relations imply the following inequalities, 



> — Ul--Ui,, (1.22a) 

g w g MM' ^ ' 

> --Ill--Ut, (1.22b) 

o MM o MM' ^ ' 



which are as strong as Weingarten's inequality. From the expressions for the S, P and T 
channels we can immediately derive the inequality 

n^>n^ + nj,^„ (1.23) 



which is stronger than Weingarten's inequahty. Because both the vector and axial mesons 
couple to the tensor current (see section 1.4) it follows that the pion mass is not only lighter 
than the scalar mesons but also lighter than the vector and axial mesons, i.e., < rrip 
and m^r < rria^ , which again confirms the Vafa-Witten theorem. 

In the chiral limit, the pion is a Goldstone mode whereas the mesons in the other 4 
channels remain massive. This requires a precise cancellation at large separations. More 
precisely, up to exponentially small corrections we have 

Qsix) = Qp{x) = \Qv{x) = \Qa{x) = ^Qt{x), for x ^ oo. (1.24) 
4 4 o 

In other words, all components Qi, in the tensor decomposition of the propagator, converge 
to the same limit at large separations. 

As a useful practical application of these general relations, we derive the corrections 
of the quark-condensate to the free propagator. In the so called 'vacuum dominance' 
approximation^, the propagator is color diagonal and simplifies to 

S = So + ^\<ijij>\. (1.25) 

If the end points are taken along the r— direction, the free propagator is given by 

^oW = (1-26) 
For the correlators in the different channels one finds 

n^(-) = ;^ + ^l<^^>l'' (i-27b) 

nj.^.(r) = -i|<^^>p. (1.27e) 



-'^The reader should be warned, that 'vacuum dominance' as introduced in |j docs not lead to the 
corrections we derive now, but rather to a kind of radiative correction to it. The reason is that we 
consider the OPE in space-time representation and include contributions that are 'regular' at a; = y, 
while in Q the OPE was studied at large momentum transfer, which includes only singular terms. See 
Oil for further discussion. 
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For completeness, let us also mention the results for the nonzero non-diagonal corre- 
lators: 

As we will see below, in all channels the sign of the corrections to the free quark 
contribution coincides with the experimental trends. Thus, although the vacuum domi- 
nance approximation does not provide a quantitative explanation of the behavior of the 
correlators it certainly gives some qualitative insight into it. 

Finally, we derive an inequality that can only be proved under a very plausible but 
not general condition. Let us consider the difference Qi — Q^. It can be written as 

Qi-Q, = \< Tr(l - 75)^Tr*(l + 75)^ >, (1.29) 

where the average is over all gauge field configurations. If the fermionic modes of opposite 
chirality are uncorrelated, this average can be factorized as 

Qi-Q, = ^ <Tr(l-75)S><Tr*(l + 75)5>, 

= ^ < Tr5 >< Ti*S >, (1.30) 

where the second equality can be established because < Tr75S' >= for a parity invariant 
vacuum state. In terms of the mesonic correlation function the positivity of Qi —Q5 leads 
to the inequality 

I^vix) > n^(x), (1.31) 
which is obeyed by our numerical results presented below. One also has the inequality 



ny(q') > n^(g), which can be proved rigorously for Euclidean momenta [23 . 

The conditions for (1.30) are realized for a random ensemble of instantons. In the 
opposite case, of an ensemble of instanton— anti-instanton molecules, we have 

TtS = Tr75S (1.32) 

before averaging over the collective coordinates. Therefore, Qi = Qs, and the axial and 
vector correlators are degenerate. 
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1.4. Spectral functions for flavored correlators 



Next we discuss the spectral functions for the different flavored correlators. They will 
be approximated by the sum of one or two meson resonances in the form of a 5 function 
and the contribution from the continuum. First, we consider the S, P-case in which the 
coupling of the current to the resonance state is deflned by the matrix element 

< 0|j^'^(x)|p >= A^'^ , ^ e-''P\ (1.33) 
The total spectral function is then given by (see ref. 
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P'^'^^s) = Xip6is - mip) + —eis - Eos,p). (1-34) 

The continuum contribution can be obtained perturbatively. However, it also follows from 
asymptotic freedom. At short distances the correlator approaches the free correlator given 
by Alternatively, the continuum contribution is given by the integral 

/ dsD{^/^,x)p{s). (1.35) 

For small x the main contribution comes from large s. In this region p{s) ~ s on dimen- 
sional grounds, and by an explicit calculation of the integral one flnds p{s) = (3/87r^)s. 
In the vector channel the spectral function is transverse and can be written as 

P:. = {-9,. + '-f)pW)- (1-36) 



Its scalar part pX^iff) is positive deflnite ||2T|. The coupling constant of the vector current 
to the p meson is deflned by 

< Q\3l{x)\pp >= i\^^e-^^- (1.37) 

where is the polarization vector of the p— meson. Using the same approximation as in 
the S", P-case the total spectral function is given by 

- pIM = ^K^i' - + |^^(^ - ^op). (1.38) 

For the axial vector channel we have the additional complication that also the longi- 
tudinal combination contributes to the spectral function. 

Pi = {-9,u + ^)p^(g^) + q,q.pi{q\ (1-39) 
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where both Pxiff) and p^iff) are positive definite. The longitudinal component of the 
current couples to the pion. The relevant matrix element is conventionally written as 

< 0|i;f (x)|7rp >= if^p^^^^e-'^^, (1.40) 

^(27r)3 

with the experimental value of /tt = 133.7 ± O.lSMeF. The spectral function in this case 
is given by 

- PtM = -flml5{s - ml) + 3A^/(. - J + —e{s - E,^,\ (1.41) 
where the coupling constant of the oi meson to the axial current is defined as 

< 0|j7(x)|aip >= iK^L=e-^^\ (1.42) 

Our last diagonal correlator is the tensor correlator. In this case the contribution 
of the free quark propagator vanishes which allows us to extract the parameters of the 
resonances in an much 'cleaner' way. The tensor structure of this current has two 3- vectors 
(like the electric and magnetic field) with opposite P— parity, and thus couples to both 
vector and axial mesons. One can introduce the corresponding coupling constants as 

<0|j^z.|p> = fpie^Pu - e,,p^) J e~^^^, (1.43a) 

< > = faien,,ai3Paef3 — e''^'', (1.43b) 

resulting in the spectral function 

pI,.(s) = 6/>J<5(s - ml) + ^iiyjis - m^J. (1.44) 

Let us now briefly mention the non- diagonal correlators. We have 5 different flavored 
currents and therefore 10 different pairs. The combinations SP, SA, PV, VA and AT 
are zero as a consequence of the parity invariance of the QCD action. The simplest way 
to show this is to take the separation between the two currents in the 4-direction. The 
desired result then follows immediately from a parity transformation in the path integral 
for the average correlator. Two combinations, ST and PT, are zero because an anti- 
symmetric tensor cannot be constructed out of only one vector. The SV correlator is zero 
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as a result of the conservation of the vector current. Formally, this can be proved from 
the Ward identities < Tuu{x)u'y^u{0) >= and < Tuu{x)d'^ ^d{{i) >= 0. 

Only the PA and VT non-diagonal correlators arc non-zero. In both cases the con- 
tribution of the free propagator to the corresponding spectral function vanishes, and 
consequently, the continuum contribution is absent. If 'pion dominance' is used we find 

P^A'i')='P^>^f^^(l'-^l)- (1-45) 
In the \^T-case, the natural approximation is to saturate the intermediate state ex- 
pansion by the p— meson. The spectral function at low momenta is then given by 

pli^Aq') = ^tpAjAi' - ^D- (1-46) 

1.5. Unflavored correlators 

We consider correlators of the unflavored currents defined in (1.3). The general struc- 
ture of the correlator of uFau and uTbu is given by 

n =< TrS{x, y)rAS{y, x)rB > + < TrV aS{x, x) TrVBSiy, y) > . (1.47) 

The characteristic feature is the appearance of a two-loop or annihilation diagram. The 
one- loop diagram is absent if the flavor content of the two currents is different. 

The vector and axial traces, Tr7^5'(a;, y) and Tr7^755'(a;, y) are not well defined for 
y — > X. Both traces have to be regularized by the inclusion of the path ordered exponential 
P exp —i A^dx/^ and the subtraction of the free propagator (which only contributes in 
the vector channel). In the axial case there will be an additional contribution due to the 
chiral anomaly. 

However, according to Zweig's rule, all flavor changing interactions arc strongly sup- 
pressed in vector and axial channels. Therefore, we will ignore the corresponding con- 
tributions in our calculations below and approximate the unflavored axial and vector 
correlators by its one-loop contributions only. 

The situation is different in the scalar, pseudoscalar and tensor channels. In this case 
the annihilation diagrams are finite. Let us consider the correlator of the current 
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with r^"^^"^ equal to 1, ^75 or a^^^ and a similarly defined current jY^^- Only the annihi- 
lation diagrams contribute: 

= - < x)r-- TrSaiv, y)r-- > . (1.49) 

According to the identity (1.16) we have 

Trr^^">S(x,x) = -Trr^"°t^(x,x), (1.50) 

which implies that the trace is real in the pseudoscalar case and purely imaginary in the 
scalar and the tensor channel. We thus find that the two-loop pseudo-scalar correlator 
is negative for x close to y whereas the contribution from this region to the scalar and 
tensor 2-loop correlator is positive. 

Next, we derive the non-diagonal spectral functions for the scalar and pseudoscalar 
fiavored currents. In the scalar case uu and dd can be written in terms of the sum 
and difference of j'^ and (see eq. (1-3)). Because the contribution of the free quark 
propagator vanishes we find the spectral function 

pS, 2-ioop(^ ^ ^ (2^f5^^q) < Q\f\Q >< Q\f\Q > +A^5(s - m^) - \]5{s - ml), 

(1.51) 

which indeed gives rise to a positive correlator at short distances (the first term in the 
r.h.s. has to be manipulated carefully). 

The ui'j^u and di'-f^d currents can be written as a linear combination of j^, and j'^ 
which allows to express the two-loop pseudoscalar correlator as 

= < 0\Tr\x)r\y)\0 > +i < 0\Tf{x)f{y)\0 > +1 < Q\Tj\x)j\y)\Q >(1.52) 

Again, the free quark contribution to the correlator vanishes which allows to write to 
corresponding spectral function as 

pP 2-ioop(^) ^ -^XloSis - mlo) + lxlS{s - ml) + ^Xl6{s - mj). (1.53) 

Let us, for simplicity, consider the limit in which all quark masses are equal. Then 
mTrO = mg and A^o = Ag. The S'C/(3)— singlet spectral function, p^' , is then given by the 
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sum 



P 




\l5{s-ml) + ^^e{s-E,p) 



(1.54) 



where the 1-loop pseudoscalar spectral function coincides with the flavored pseudoscalar 
spectral function obtained in eq. (1.34). Several important conclusions can be drawn 
from this result. First, the continuum threshold in the r/'— channel is the same as for 
other pseudoscalar mesons. Second, the r/'— mass is determined by the two-loop spectral 
function, whereas the mass of the other pseudoscalar mesons follows form the 1-loop 
spectral function. In the r^'— channel the contributions of all other pseudoscalar mesons 
cancel because of a detailed compensation between 1-loop and 2-loop contribution. At 
relevant distances of the order of 1 fm the pion correlator is several orders of magnitudes 
larger than r^'— correlator, and therefore both 1-loop and 2-loop contributions have to be 
calculated extremely accurately in order to determine the mass of the rj' particle. 

Finally, we show that in the quenched approximation the positivity of the 77'— correlator 
is violated. The simplest argument goes as follows. Suppose, that all quark masses are 
equal. Then the singlet 77'— correlator is given by the sum of the one-loop pseudoscalar 
correlator and Nf times the two-loop correlator 11^^^. In the quenched approximation both 
contributions do not depend on Nf, and it is clear that if the two contributions balance 
each other at one value of A^/, the correlator becomes negative for Nf + 1. 



Before going into the numerical calculations of multi-instanton effects, let us briefly 
remind the reader of the qualitative effects of the one-instanton gauge field fluctuations. 
The main physical phenomenon, tunneling through a barrier that separates gauge fields 
of different topology, is related to a rearrangement of the light quark states: some of them 
'dive into the Dirac sea' during this process, whereas others 'emerge' from it. This process 
is described by the so called 't Hooft effective Lagmngian 



1.6. Qualitative role of the instanton-induced interactions 



LeS ~ n(?/^o)(^Og/) 
/ 



(1.55) 



16 



where the quark fields qf of various fiavors (/ = d, s) are projected onto the so called 
fermionic zero mode ipo{x), which is the solution of the Dirac equation DiIjq{x) = in 
the field of an instanton. The zero modes play the role of wave functions of the states, in 
which the quarks are produced or absorbed during tunneling. They depend in a known 
way on collective coordinates of the instanton, its position, size and color orientation. 

An important observation is that the chirality of the zero modes is directly related to 
the topological charge of the gauge field: there is only a left-handed fermion zero mode for 
an instanton and only a right-handed fermion zero mode for an anti-instanton. For anti- 
quarks, the opposite holds true. As a result, if one ignores quark masses, it is impossible 
to close the loop and return to the emitted quark because its chirality flipped on the 
way. Therefore, in the chiral limit, 'solitary' instantons are absent. They only appear in 
groups of zero total topological charge, exchanging produced quarks among themselves. 
Obviously, in order to break chiral symmetry, such 'group' has to be infinitely large, 
percolating the entire space-time volume of the Universe. The low-lying fermion modes 
become 'collectivized', and occupy the entire space-time as well, instead of being localized 
near a particular instanton. 

At this point, one can ask the question: in what sense the 'one instanton approxi- 
mation' can be formulated, if instantons cannot be isolated and the low-lying fermion 
modes are delocalized? Well, provided that the instantons still form a reasonably dilute 
ensemble, the delocalized modes can be written as a superposition of localized zero modes, 
■0a(2^) — J2n Cx'^o{x — Zn), whcrc Zn is the center of the instanton. The coefficients behave 
as C" ~ l/V^, but, because '4>{x) ~ 1/x^ for large x, only nearby instantons contribute 
to |'(/'A(a^)P, and for a reasonably dilute ensemble |-?/'A(ic)p is peaked at the center of the 
instantons. For the same reason, provided the end points x and y (for which correlator is 
evaluated) are located within the radius of a single instanton j, the quark propagator is 
dominated by the term [^Pq{x — Zj)ipl{y — Zj)][J2x \Ci\'^/{X + im)]. The first factor, pro- 
jected on quark fields, reproduces the original 't Hooft interaction. At the same time, the 
second factor, (which for a single instanton in the perturbative vacuum is 1/m (because 
there is only one zero mode) is now different: it has a non-singular chiral limit resulting 
from the many different modes that contribute to the sum over A. Since this factor does 
not depend on x — y, it and can be ignored for qualitative studies of instanton effects, as 
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if only one instanton was present. 

A phenomenological discussion of the instanton-induced four-fermion interaction was 
first made for the pseudoscalar channels in 0, and similar studies have eventually led 
to the 'instanton liquid model' [|10|, [0]. Below we make it even more convincing, by 
considering instanton-induced effects in the scalar channels as well [Q. 

The specific flavor and chirality-flipping structure of the 't Hooft interaction leads to 
the following important conclusions: 

1. The effect of the interaction is of first order for the scalar and pseudoscalar correlators, 
hut is absent in the vector and axial channels. 

2. The sign of the corrections is opposite for scalar and pseudoscalar channels. 

3. Due to its flavor structure, uudd, the one-instanton corrections have opposite signs for 
isospin 1 and correlators, which are essentially given by {uu ± dd)"^ . 

The first point agrees with phenomenological observations ||12|, according to which 
much stronger deviations from the asymptotic freedom are seen in spin-zero channels 
than in spin-one channels. The last two statements provide a qualitative understanding 
of the short distance behaviour of all spin-zero correlators. Ignoring for simplicity the 
strange quarks, we have four channels with different parity and isospin, which we denote 
as Tx (/^ = 1'), i (/^ = 0-), a (/^ = 0+) and 5 (/^ = One observes that in the 
one-instanton approximation the corrections in the vr and a channels have the same sign 
as free propagation, which implies a larger correlator and an attractive interaction. The 
opposite, with a repulsive interaction, is found in the two other cases, the rj' and the 6. 
This behavior of these four channels is indeed seen in the real world, which shows that, 
by studying instanton-induced interactions, we are on the right track. 

In principle, one could proceed along these lines, in which the 't Hooft interaction is 
treated as an effective multi-fermion interaction that can be taken into account system- 
atically in perturbation theory. However, we do not follow this approach, because we 
will include all such diagrams by an exact diagonalization of the Dirac operator in the 
subspace of zero-modes. 



2. The 'flavored' mesonic channels 
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In this and the following section we present our numerical results obtained for an 
ensemble of 128 instantons and 128 anti-instantons in a periodic box of (3.36)^ x 6.72 fm^. 
The end points of the correlator are taken along the longest axis, in which way most of 
the finite size effects due to periodicity are suppressed even for the maximum considered 
distance of 1.615 fm. The average correlators have been obtained by averaging over 50 
different gauge field configurations, distributed according to the invariant measure of the 
integration over the collective coordinates, and, for each configuration, by averaging over 
100 different points x for a fixed value of x — y. The size of the instantons is taken fixed 
at p = 0.35/m. 

The definition of the quark propagator in the 'random instanton vacuum' was ex- 
plained in detail in Let us only comment here that for the present calculations the 
quark masses have the value of 10 MeV for u and d, and 140 MeV for s. The former is 
still larger than the physical mass0, but it is about as small as we can tolerate in order to 
be insensitive the the artifacts in the spectrum of the Dirac operator at small virtual- 
ities. The latter is not fitted to any particular quantity but is taken as a 'common sense 
value'. 

In order to check the volume dependence of the correlators, we also have performed 
calculations with fewer instantons. A significant volume dependence is only found for the 
pion and, to a much lesser extent, for the kaon. For example, decreasing the number of 
instantons from 256 to 64 at constant density, raises the pseudoscalar correlator at 1.5 
fm by about a factor of 2.5, whereas 2/3 of this factor occurs by decreasing the number 
of instantons from 128 to 64. Under the same conditions, the kaon correlator changes 
by only 25 percent. The volume variation from decreasing the number of instantons by 
a factor of two is inside the error bars for all other mesonic correlators that have been 
studied. 

2.1. Two pseudoscalar channels: n and K 



^ Therefore, for a physical quark mass of rriu + rud = 11 MeV , our pion mass should be m^j' 

experimental /■2Q 
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In this section we consider correlators of the flavored currents 

jn = -^{ui'j^u- di-f^d), (2.1a) 
jK = ui-i^s. (2.1b) 

It is obvious, that as these pseudoscalars are the lowest excitations of the QCD vacuum, 
the Goldstone modes, they can tell us a great about its long-range structure, like for 
example the chiral condensate. However, unbeknownst is that these channels provide 
very significant information about the short-range vacuum structure as well. 

Another, but closely related, exceptional feature of the pseudoscalar correlators is 
that all components of the propagator enter with the same positive sign (see eq. (1.18)). 
Therefore, they are the easiest to measure: no cancellations take place. For the same 
reason, though, the systematic error due to the finite volume is largest. 

In Fig. 1, our results for the vr and K correlators are compared with phenomenological 
expectations. As usual, they are normalized with respect to the free quark correlator, 

no(x) = (2.2) 

which corresponds to the simple 1-loop diagram. Therefore, the fact that our points 
approach 1 at small distances is a general consequence of 'asymptotic freedom', and 
should be independent of the non-perturbative vacuum fields. What is studied is the 
deviations from this simple behaviour. 

The solid (phenomenological) curves have been obtained from the expression 

pfl . A^B,™,, + r ,EE^o(E, .)(1 + (2.3) 

IIo(x) 3 4 JSo TT 

that can be obtained from the spectral function (1.34) and the free propagator (see 
for an explanation of the perturbative correction). The parameters involved, namely the 
coupling of the pseudoscalar current to these particles, and the 'continuum threshold' Eq 
above which the quarks are assumed to be produced as free ones are given by: 

= (480Me1/)^ A^^ = 1.24A^, Eo=lAGeV. (2.4) 

Unfortunately, the uncertainty in the absolute value of quark condensate translates into 
an about 50% uncertainty of the couplings. The uncertainty in Eq is important only in a 
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very limited range of distances of about 0.3 — 0.5 fm, where resonance and non-resonance 
contributions are comparable. For a more detailed analysis of these parameters we refer 
toref. g. 

The dashed line, shown for comparison, represents the 'vacuum dominance' approx- 
imation (see above). Although this correction has 'good intentions', the correct sign, it 
describes neither our points nor the experimental curve^. More elaborate OPE formulae 
fail in these cases as well (see [0, and 

The parameters of (2.3) without the perturbative correction can also be obtained from 
a fit to our data points. The results are shown in Table 1. For the pion both the mass 
and the coupling constant agree perfectly with experiment. For the K meson the mass 
is right where it should be, but the coupling constant is too small. For a comparison to 
other channels and a general discussion of the parameters we refer to section 4.2. 



2.2. Vector channels: p, K* and </> 

We use vector currents with the same normalization as in []T|Q 



= -^[wy^u - d'j^d] or wjf.d, (2.5a) 

= -^[uj^u + d^^d], (2.5b) 
Jt = sj,s, (2.5c) 
j^* = uj,s. (2.5d) 

The phenomenological expectations of the correlators can be obtained from the spectral 
function (see section 1.2). However, in order to make contact with experiment, one uses 
that the spectral density is proportional to the ratio 

Ri{s) = (Te+e-^i{s)/ae+e-^^L+f,-{s) (2.6) 



■^Note, that in this case the disagreement with experiment is not affected by the uncertainty in the vaiue 
of the quark condensate, because both the contribution of the pion and this correction are proportional 
to its square. 

^Warning: other authors, e.g. 0, use a different normalization. Any definition that includes the 
charge of the quarks is confusing if one compares different channels. 
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of the e~^e~ annihilation cross section into hadrons with quantum numbers i and the 
cross section for muon pair production. Neglecting the muon mass, the latter is just 
ae+e-^^+fi- = (47ra^/3s). The trace if the spectral function in the vector channel is 
related to this ratio by 

2npUs) = (2.7) 

which after substitution in the dispersion relation leads to 

< 0|Tj;(x)j;t(0)|0 >= / dsDis'/\x)sR''{s). (2.8) 

This relation can be used to extract the 'experimental definition' of the correlation func- 
tion. 

Information on strange vector channel, marked by i^'*(892), is derived from the weak 
decay process r — > z/^- + hadrons due to Cabbibo mixing. We also include the double- 
strange (j) channel. Strictly speaking, it is related to both one-loop and two-loop diagrams, 
but according to experimental observations ( the famous 'Zweig rule'), all flavor- changing 
transitions are strongly suppressed in vector channels, and the two-loop contribution 
cannot be a significant correction. All final states with a kaon pair are included (see 
details in [|^]) in the phenomenological curves. Also the u and p correlators only differ 
by strongly suppressed flavor changing two-loop diagrams. Within our accuracy this 
difference is not visible, and therefore we only show results for the p— correlator. 

Numerical results for the vector mesons are shown by the solid lines in Fig. 2. As usual, 
correlators are plotted in a normalized way, namely as n^^(x)/n°^(a;), where 11°^ (x) = 
— G/vr^x^ corresponds to the simple loop diagram, describing free propagation of a pair of 
massless quarks. 

The general striking observation in all vector channels is that the contributions 
to the spectral function of the lowest meson and other states complement each other in 
such a way, that the ratio n(x)/n°(x) remains close to one up to the distances as large 
as 1.5 fm\. This phenomenon, called superduality, generalizes a well-known 'quark-meson 
duality' argument to much larger distances. From an experimental point of view it means 
some 'fine tuning' of the parameters of all vector states. On the other hand, in the 
language of field theory it implies that in the vector channels, for some (so far unknown!) 
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reason^, all non-perturbative corrections to the free quark propagators miraculously cancel 
each other, until the correlator has dropped by a few orders of magnitude. 

Results of our calculations are compared with the phenomenological expectations (solid 
curves in Fig. 2) obtained from the experimental i?— ratio as discussed aboveQ First of 
all, 'superduality ' is qualitatively reproduced by the calculated points. Keeping in mind 
that the correlators drop by several orders of magnitude over the range of distances under 
study, even the quantitatively agreement is actually very good. The fact that all our 
curves are somewhat below the experimental ones came as no surprise: in the 'instanton 
vacuum' under consideration there are no radiative corrections to correlators. To first 
order, they are known to be (1 + as{x)/TT + ...), which means an increase in the value of 
the correlator value by about 10 percent. If the next terms produce a smaller or similar 
correction, the agreement with data is actually nearly perfect. 

Let us also comment that the splitting between the three vector channels of flavor 
content ud, us and ss is surprisingly small up to x = 1.5 fm. As the strange quark mass 
is generally far from being negligible {e.g. in the way it affects the propagators), this 
observation means that the 0{ms) terms have the tendency to cancel among themselves. 
The theoretical reasons for this are unclear. At large distances, x ~ 1.2 fm, a splitting 
between the different vector mesons appears, but the order is wrong: in the instanton 
vacuum the p meson is heavier than the meson. One may also look at our calculated 
results in a different way: instead of comparing them to the 'phenomenological correlator', 
they can be fitted directly to some parametrization, which allows us to extract some 
'hadronic parameters'. As in the previous section, we use the standard three-parameter 
expression, which now looks as follows: 



The coupling constant is defined as in eq. (1.37) which differs from the definition that 



comparison between the various channels can be made most unbiasedly. 

^ As was discussed briefly in section 1.6, a partial explanation of this phenomenon is provided by the 
instanton-based theory: at the one-instanton level there is no correction to vector (and axial) correlators, 
as opposed to scalar and pseudoscalar ones. 

^The 1-loop radiative corrections are included for the p and the ai channel, but are absent in case of 
the (f) meson, see [|]. 

^In 10] dimensionless coupling constants are used, which are related to ours e.g. as Xp — ^/^rri^/g^ 




(2.9) 



is used for extracting coupling constants from the dataQ. We think that in this way the 
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In particular, in |^ we have used the matrix elements of electromagnetic current 



< Oli;:-'^- Imeson >= f^^^m^esone^ (2.10) 

which then lead to a universal relation containing a well measured quantity, the electro- 
magnetic branching width: 



Smmcsnnr (meson e+e" 



/meson ^^o ■ { ■ ) 



As the electromagnetic current contains factors related to the electric charges of quarks 

we have to remove them which, in absence of mixing between the different channels, leads 
to the coupling constants 



= V2fpmp= {A09±5MeV)^, (2.2a) 
A^ = 3V2f^m^ = (390 ± 5 MeVf, (2.2b) 
= 3f^m^ = {A92±15MeVf. (2.2c) 



The K* coupling constant, fx*, is defined in the same way as fp and can be obtained 
from the weak r— decay (see fll). For X^* one finds 

.rriK* = (448 ± 25MeV)'^. (2.14) 

In table 1 these accurately known coupling constants are compared with our calculated 
results. 



2.3. Axial correlators: ai, Ki and /i 

Currents and correlators in these channels are defined in exactly the same way as 
for the three vector channels in the previous section with the obvious substitution 7^ — »• 
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j7 = ■^[u^f.lbU - d7^,l5d] or U'jf.-f^d, (2.15a) 
= s7m75S, (2.15b) 
= uj,j,s. (2.15c) 

Again, we assume that Zweig's rule is very accurate so that annihilation diagrams can 
be ignored in the ss-type channel. 

However, as we have seen in section 1.2, there is a significant difference: the axial 
correlator correlator 11^^ (g) has a non-zero longitudinal part, proportional to q^qu- For 
zero quark masses it only contains the contribution of massless Goldstone modes, pions 
etc.. For non-zero quark masses, axial currents are not conserved, and a massive pion 
contributes to the spectral function. Thus, generally speaking, in the axial cases one 
has two different correlation functions which can be split in, for example, transverse and 
longitudinal parts. 

In this work we concentrate on one particular combination: the trace H^^. The reason 
is that at this moment we are mainly interested in the ai meson, and not in the pion 
which, in this combination, is suppressed by the small parameter and is unimportant 
at not too large distances. The details of the extraction of the relevant ai parameters from 
the T lepton decay data can be found in [Q]. The normalized correlator can be represented 
as follows^: 

= 3XlD{m^,,x)^ - />^D(m.,x)^ + ^ J^^ dEE'D{E,x), (2.16) 

where the last term represents the non-resonance contribution. It is given approximately 
by its asymptotic form with some threshold. Unfortunately, the r lepton is not heavy 
enough to allow its determination, and we use two 'reasonable' values, -Eo=l-5 GeV and 
Eq =1.7 GeV, to show the sensitivity of the correlator to this parameter (see two solid 
curves in Fig. 3). Also in this case the 1-loop radiative corrections (not displayed in 
(2.16)) are included in the phenomenological curves. 

®The reader should be warned, that the sign of the pion term in this formula in is wrong. As a 
result, this correlation function changes sign and becomes negative at large distances, where the pion 
term becomes dominant. 
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Our calculated results for the three flavored axial vector channels are presented by 
the data points in Fig. 3. As for the channels considered above, the agreement between 
calculations and phenomenology is amazingly good. All general features, including the 
fact that axial mesons are significantly heavier than the corresponding vector mesons, are 
obviously reproduced. Moreover, the numerical agreement with expectations at small x 
would be nearly perfect, if radiative corrections would have been taken into account. 

Let us also comment, that we find surprisingly close results for ud, us and ss type 
correlators in a wide range of distances, until approximately x = 1 fm. It means that, as 
in the case of the vector channels, in the whole region the 0{ms) terms have a tendency 
to cancel among themselves. 

At large distances, above 1 /m, the contribution of the much lighter pseudoscalars 
(which we tried to suppress, by taking the trace) still becomes dominant. Somewhat 
surprisingly, the data look good enough for the extraction of their coupling constants and 
masses. Using the expression given above we have made a five parameter fit with results 
shown by the dashed-dotted lines in fig. 3. 

The fitted parameters of axial mesons and non-resonance continuum are listed in Table 
1 and are discussed, together with others, in section 4.2. Here we only present the results 
and some comments concerning the pseudoscalars. 

For the pion the fitted values, corresponding to minimal (and to the curve in Fig. 
3), are as follows: 

= 252 ± IbMeV, U = llOMeV. (2.17) 

The latter value can be directly compared with the well known experimental result of 131 
MeV. However, in order to compare the calculated pion mass to its experimental value, 
one should make an adjustment for the fact that we cannot take as small quark masses 
as required by the real world. We remind that our value for (m^ + nid) caic—'^O MeV , 
while experimentally (m^ + nidjexp ~ H MeV, and therefore we should extrapolate the 
calculated pion mass to 222^11/20 MeV = 178 ± 10 MeV which is not far from the 
observed value in the pseudoscalar channel of 142 ± 14 MeV. In other words, the model 
does reproduce the right magnitude of the mass- independent coefficient m^/(m„ -|- m^) 
for small enough quark masses. 
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The corresponding fitted values of tlie parameters in (2.16) for tlie kaon are: 

niK = 464 ± 15MeV, /k = lOOMeV. (2.18) 

Tlie mass is right on the experimental value, but the decay constant is about a factor 2 
smaller than needed. 

For the 'purely strange' pseudoscalar rjs we find from the fit: 

m^^ = 553 ± 15MeV, ^ = 90MeV. (2.19) 

As this combination is actually some mixture of rj and t]', with the two-loop effects so far 
ignored (see below), we do not think these results have a direct physical meaning, but we 
present them for completeness. 



2.4. Two 'flavored' scalars 

At this point we would like to make some remarks on the accuracy of our calcula- 
tions. All calculations have been performed in the quenched approximation where the 
small eigenvalues are not suppressed by the fermion determinant. In the one-instanton 
approximation the vector and axial correlators have nevertheless a finite chiral limit, 
and we expect that the enhancement in the spectrum for small virtualities does not 
significantly alter the correlator. However, in the case of the scalar, pseudoscalar and 
tensor correlators, the chiral limit cannot be taken in the one-instanton approximation, 
and therefore we should not be surprised if we are unable to reproduce the experimental 
results with the physical values of the light quark masses. Because the pseudoscalar corre- 
lator is the sum of all 'hidden components' the relative error is still small, and, in section 
2.1, we have found that the phenomenological curves can be obtained with somewhat 
larger quark masses. On the other hand, in the case of the flavored scalar correlators, 
the first-order instanton-induced interaction is 'repulsive at small distances (see section 
1.6). At larger distances, other corrections {e.g. multi-instanton ones) become important, 
but these correlation functions are small as compared to those considered previously. At 
a distance of the order of 1 fm a cancellation of a couple of orders of magnitude has to 
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take place between the spin flipping and the spin non-flipping components of the propa- 
gator. We necessary will find a large systematic error not only because of the anomaly 
in the spectrum at small virtualities, but also because wc use approximate propagators. 
In the tensor channel the instanton induced interactions are attractive. Consequently, 
we have much less cancellations and much more accurate results can be obtained. The 
non-diagonal correlators are less singular in the chiral limit (see, for example, the vacuum 
dominance results) , and also in this case we expect smaller errors than in the scalar case. 

After these warnings, we are ready to confront our results for the scalar correlator. 
Our results for the ud and us channels are shown by the data points in Fig. 4. One 
observes that our data decrease more rapidly than in any other channel considered above. 
In fact, our data even 'overshoot' this trend. The 'random instanton vacuum' produces 
a too strong repulsion in these channels: the correlators become negative (since we do 
not trust the results for the correlator in this region we did not even plot those points). 
This violates the positivity of the scalar correlator and it points to a defect of the model. 
Indeed, as we have stressed above, all calculations have been performed in the quenched 
approximation. 

Apart from the ao(980) particle (which is believed to be mainly an rjn or a KK system 
with flavor composition essentially given by udss), no flavored scalar mesons are listed 
by Particle Data Group, so, probably, they do not exist at all! However, continuum 
(multi-meson) states should of course exist, and, for comparison, we have plotted some 
'expectations' in Fig. 4. We have assumed the absence of resonances and have pushed 
continuum threshold up to = 2 GeV, which we consider as a kind of its upper limit. 
Therefore, the solid curve represents the most rapidly decreasing correlator, on the edge 
of what seems reasonable. 



2.5. Tensor and non-diagonal correlators 

This section deals with a set of correlators which are only rarely discussed in literature, 
e.g. they are not mentioned in reviews on the subject. However, as wc arc going to show 
shortly, they are actually quite useful, because all of them can be described very well by 
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the dominance of one intermediate state. Therefore, they are much more suitable for the 
determination of hadronic masses than all standard correlators. 

The first example is the correlator of two 'tensor currents' j^i, = ua^j^d (and also 
its strange analog, d ^ s). As the two indices are anti-symmetric in this case, the 
intermediate physical states are not the usual tensor mesons, of course, but rather some 
specific polarization states of vector and axial mesons. 

Contrary to all other diagonal correlators considered above, these tensor correlators 
do not have a free-quark contribution at small distances: the corresponding loop diagram 
vanishes kinematically because of the Dirac trace alone. (Therefore, one hopes to obtain 
the contribution of resonances in a much 'cleaner' way.) The coupling constants were 
defined in (1.43), and, in this case the normalized correlator simply follows from eq. 
(1.44) with result 

= 2vrV/XD(m„x) + 2vrV/f D(m,„x). (2.20) 

iio {X} 

Our data are shown in Fig. 5 together with an approximate fit, corresponding to 
fp = 100 MeV, fa^ = O.OMeV and experimental p and ai masses. Thus, this correlator 
is completely dominated by a single rho meson pole. 

Our next topic is the non-diagonal vector-tensor correlator. This is the only non- 
diagonal correlator that was considered in the literature in connection with QCD sum 



rules [|T^, |14|, |T5| . It is related to the baryon magnetic moments and other electromagnetic 
properties of hadrons. In particular, the susceptibility of the quark condensate xIq"^) is 
defined by the small-g limit of the integral 

J d^xe'^^U''/ = Stxiq') < qq > q,. (2.21) 

The correlator can also be expressed in terms of the spectral function. When substitute 
the p— meson dominance formula (see eq. (1.46)) in eq. (1.9) we find the sum-rule 

< qq > [q + rrip) 

For phenomenological values of rrip, Xp and Xp and our fitted value of fp we find x(0) = 
—2.0GeV^^ which differs from the corresponding value x(0) = —3.3GeV^^ obtained in 
T^ . This deviation originates in the coupling constant fp which was estimated with the 
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help of vacuum dominance in [|14|, i.e., \pfp = 2 < qq >, resulting in fp = 165 MeV 
as opposed to our value of 100 MeV. A more elaborate QCD sum rule calculation 
taking into account two additional p— meson resonances gives a slightly larger value 
of Xe.m. = —^.AGeV'"^ (at the normalization point /io = 1 GeV). 
More exphcitly, they used the spectral function 

p^(g^) = SiqplCpSiq'^ - mj) + Cp/(5(g^ - mj) + Cp>>S{q'^ - nip,,'^)], (2.23) 

where the three terms contain the contributions from the three lowest p resonances, with 
masses of 780, 1250 and 1600 MeV, respectively. The constants are related to the coupling 
constants and are equal to Cp = 2.9| < qq > \, Cpi = — 1.1| < > | and Cp// = 0.6| < qq > \. 
The ratio of the corresponding VT correlator to the free pseudoscalar correlator 11^ given 
by (for xq > 0) 

npj^(x) ^ ^4^6^^ dD{mp,x) ^ ^ dD{mp>,x) ^ ^ dD{mp>,,x) ^ 
XIq (3^) doCy doCy doCy 

is represented by the solid curve in Fig. 6a. Let us remind, there are no free parameters 
and that we have an absolute normalization. In the same figure, these predictions are 
compared with our results. One should note that although the two calculations are 
completely different in nature the agreement is quite reasonable. 

At small distances the agreement is exceptionally good. Of course, this is not acci- 
dental: QCD sum rules are based on the OPE, and should do better in this region. In 
fact, the OPE dictates that the short distance behaviour of 11^-^ is given by the vacuum 
dominance formula (1.28b). Therefore, the excellent agreement at very small distances is 
essentially do to the fact that we reproduce the right magnitude of the quark condensate. 

A natural generalization of electromagnetic susceptibility of the quark condensate to 
the case of a weak external axial current \ T6\ leads to another non-diagonal correlator: the 



pseudoscalar-axial correlator. It can be used to evaluate the axial coupling constants of 
various baryons (see also '^^), but presently we are interested in the correlator by itself. 



One can certainly argue that in this case one knows at least its long-range part, which 
follows from the pion intermediate state. Moreover, there are no unknown constants 
because the coupling of the pion to axial current is nothing else but the famous pion 
decay constant (see eq. (1.38)). Using 'pion dominance' for this non-diagonal correlator 

30 



we find the expression 



-Kfn 0{xo) + Kfn ^(-2^0) 

IJ, fj, 



(2.25) 



U^{x) 3 

In Fig. 6b, we show the numerical data together with the above formula with a fitted 
value of which shows agreement at all distances. If one takes = 191 MeV (bcause 
TTT'u + TTT'd — "20 MeV and Ayr = (blOMeV)"^, as fitted from pseudoscalar correlator, the 
result for pion decay constant reads 

f^^llQ±17MeV (2.26) 

which is the same value as obtained from the axial channel. 

The final remarks of this chapter deal with scalar-vector correlator. Our calculations 
show results that are consistent with zero at all distances. The difference between SV 
and PA correlators is seen from the following argument. Suppose that some intermediate 
scalar (or pseudoscalar) particle contributes to it, then we have a contribution of the 
type (2.25). Let us now take another derivative 9^ of the correlator: when it acts on 
the propagator it produces the mass squared of the intermediate state. But, in the chiral 
limit, both the vector and axial currents are conserved: the result should be zero. The 
dilemma is then as follows: (i) the coupling constant is zero, or (ii) the particle mass is 
zero. In the PA channel the second alternative is reahzed, whereas in the SV case only 
the first one is possible. 

Therefore, a zero SV correlator implies that the vector current is conserved in our 
calculations (see section (1.4)). Of course, this holds true rigorously if a complete set of 
states is used in the propagator. In practice, we work with an incomplete set of states, so 
the vanishing results are actually a non-trivial test of our calculations. 



3. Correlators of some 'unflavored' currents 
3.1. Pseudoscalar channels rj and rj' 

As explained above, both one-loop diagrams and 'two-loop' (or annihilation) diagrams 
should be included in the 'unflavored' channels (see section 1.5). It implies that the 
propagators should be defined at zero separations between its end points, or that the 
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regular part of the propagator should be separated from the singular part. Generally 
speaking, this is non-trivial: for example, the regular part of the axial current should 
reproduce correctly the axial anomaly. Fortunately, significant simplifications take place 
for pseudoscalar currents because the Dirac trace kills all chirality non-fiipping terms in 
the propagator. In practice, the 'double-loop' diagram was calculated from the zero-mode 
part of the propagator only, which is non-singular by construction. 

Let us first explain the definition of the correlators. For two massless fiavors, we can 
distinguish two different components 

A{x) = < Tuij5d{x)dij5u{0) >, (3.1a) 

B{x) = < Tui'y5u{x)di'j5d{0) >, (3.1b) 



related to one and two-loop diagrams respectively. Inclusion of a massive s quark brings 
in three more functions, 

C{x) = < Tui^5u{x)si'y5s{0) >, (3.2a) 
D{x) = < Tsi-f5s{x)si-f5s{0) >, (3.2b) 
E{x) = < Tui-f5s{x)si-f5u{0) >, (3.2c) 



with D{x) possessing both contributions, and C{x), E{x) only one of them. All 5 functions 
can be obtained from correlators in the four light mesonic channels tt, t] and rj', plus 
one non-diagonal correlator r) — rj'. 

Although we do not really use the strange quark mass as a small parameter, we still 
use the SU{3) singlet and octet currents defined in eq. (1.3). In this section we focus on 
the following 3 correlators 

< Tj,{x)n{Q) >= ^{2A + AB + AC + D), (3.3a) 

< Tjsix)jsiO) >= ^{A + 2B-4C + 2D), (3.3b) 

- < Tj^{x)jsiO) >^^{-A-2B + C + D). (3.3c) 



W{x) = 
n««(x) = 

U^\x) = 
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Note, that in the SU{3) symmetric case (m^ = = m^), B = C and D = A + B. 
Therefore 11®^ = A and 11^^ = 0, as they should: the t] and vr"*" channels are the same, and 
there is no mixing. However, due to the chiral anomaly, 11^^ = A + 3B ^ U^^. 

In Fig. 7, our results are compared with the phenomenological expectations The 
diagonal 1] and rj' correlators are shown in Fig. 7a. The obvious observation is that while 
the model is doing an excellent job in reproducing the i] curve, it certainly 'overshoots' 
in the t]' case: the positivity condition of the correlator is violated. In the rj case the two- 
loop contributions nearly cancel among themselves (exact in the chiral limit), whereas all 
one-loop contributions add up with the same sign. The 77'— meson remains massive in the 
chiral limit which requires a delicate balance between the positive one-loop and negative 
two-loop contributions (see section 1.5). The latter are closely related to the correlations 
of the topological charge density, which are not implemented properly in the 'random 
model'. 

In Fig. 7b we show our results for rj — rj' mixed correlator. Let us remind that idea 
of mixing was introduced in the days of non-relativistic quark model. One can imagine 
that the singlet and octet states rji and ?78, produced by the singlet and octet currents at 
low energies, are a mixture of the physical rj and rj' states. The non-diagonal correlator 
should look as 

= XrjXri' cos 9 sin 9[D{mr^, x) — D(m^/,a;)] + (mixing in the continuum). (3.4) 

The curve shown for comparison uses A,, = 1.3A,r, A^/ = 0.7A7r and a angle of = 20°. 
The sign and qualitative behaviour is right, but the angle (or couplings?) seems to be 
larger than expected. However, much more work is needed in order to understand this 
mixing phenomenon properly. 



This channel is different from all others considered above, because the corresponding 
current 



3.2. The isoscalar scalar (or a) channel 




1 



(3.5) 
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has a non-zero vacuum expectation value. As a result, the correlator in question has a 
non-zero disconnected part 

W{x) ^2\<qq>\^ (3.6) 

that dominates at large distances. 

Therefore, our data are plotted in two different ways. In Fig. 8a we show the mea- 
sured signal, normalized, as usual, to the perturbative correlator. The solid line shows 
the phenomenologically expected contribution of 2tt^\ < ipip > px^/3 (for < ipilj >= 
— (240MeV")^). One observes that the correlator exceeds the perturbative correlator by 
a big factor, and that our points are somewhat above the phenomenologically expected 
(solid) curve. The reason is that our model produces a slightly larger value of the con- 
densate, namely, | < ipip > \ = 2.22/m~^ = {257MeVY, which perfectly agrees with the 
large distance behavior of our points. 

However, the disconnected signal is not very interesting by itself. The physical excited 
states are related to the 'connected' part of the correlator. Unfortunately, as shown in 
Fig. 8b, the accuracy at large distances is greatly reduced after the subtraction of the 
disconnected part. It is important to note, that the connected part of the correlator is 
still quite large. This means that the effective interquark interactions are indeed strongly 
attractive. Of course, this is hardly surprising because in the one-instanton approximation 
this should be the case (see section 1.6). 

As can be seen from the dashed lines in Fig. 8b, the connected part can be very well 
fitted by our standard three-parameter fit, resonance plus continuum as in (1.34). The 
following values of the parameters were obtained: 

m<, = 543MeV, K = {bOOMeVY, = 1160 Me\/ (3.7) 

Those parameters show that the 'instanton vacuum' contains the 'sigma mesonf^ the 
famous enhancement in the tttt cross section which plays a prominent role in nuclear 
physics and is the basis of the 'sigma model' of Gell-Mann and Levy. 

Moreover, the particular parameter determined best of all is the sigma coupling con- 
stant. It is interesting to note, that it happens to be the largest of all mesons. It implies 

^ Of course, on the basis of these data one cannot answer the old question, whether the 'sigma-meson' 
is a resonance or just a well-correlated tttt pair, but one definitely sees that in order to reproduce our 
results for the correlator one should have a peak in spectral density around 600 MeV. 
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a very compact sigma particle, obviously due to the very strong attractive interaction in 
this channel. 



4. Discussion 
4.1. Comparison with other works 

In the discussion above we have compared our results for correlation functions with 
phenomenology. Now, we are going to compare them with other theoretical results. There 
generally are three different categories: (i) the operator product expansion (OPE) in the 
context of QCD sum rules with many parameters (or 'condensates'); (ii) other instanton- 
based calculations; and (iii) lattice QCD calculations. 

We will not discuss the OPE-based expressions, partly because it was done in detail in 
and partly because we already discussed its simplest version (based on the quark con- 
densate contribution) in section 1.4. For the tensor and non-diagonal correlators we have 
given a more elaborate discussion because no other sources of information are available. 

The instanton-induced effects were originally discussed in |10[ to the first-order in 
instanton density. In section 1.6 we have given a qualitative discussion which is even 
somewhat wider than those works, because it also includes the scalar channels. 

Further developments include the analytic attempts to account for the multi-instanton 
effects by means of summing a subset of diagrams |TB], The first numerical 



study of correlators in the interacting instanton vacuum p6| showed surprisingly good 
results. We do not discuss them here for two reasons, (i) Within their much larger 
statistical and systematical uncertainty the old calculations agree with the present ones, 
(ii) A much more accurate and detailed study of the 'interacting instantons' vacuum in 
under way and will be discussed elsewhere [|19| . 



In the last part of this section we compare our results with the first lattice calculations 
of point-to-point coiielaXion functions, which were recently reported in Results for four 
non-strange channels, shown in Fig. 9, show a qualitative agreement in all cases. In the 
vr and p channels one even finds, somewhat curiously, that the phenomenological curve is 
in between the two calculations. The axial channel agrees very well, while in the isovector 
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scalar (S) case our results definitely 'overshoot' the repulsion (see discussion above). 

As one compares these two sets of data, resulting from completely different theoretical 
approaches, one should keep in mind that the plotted curves do not really represent the 
measured signal: for convenience of representation the correlators are normalized with 
respect to the perturbative correlators. Actually the measured correlators change by 
several orders of magnitude, so their deviation at the largest distances by, say, a factor of 
two is in fact a surprisingly good agreement. As we will show in our next paper of this 
series this agreement is observed in the baryonic channels as well. 

The reason for the intriguing agreement between the two set of correlators can and 
should be investigated further. The most obvious way to proceed is, of course, to 'cool' 
the lattice configurations until one is left with mostly instantons, and then to redo the 
measurements of the correlators. If the resulting correlators happen to be about the same, 
our main point the dominant role of the instanton-based quark interactions in the QCD 
vacuum, would finally be demonstrated explicitly. 



4.2. The fitted masses and coupling constants 

In our discussion above we have concentrated on correlation functions as such and 
compared results directly with phenomenological curves [|I| based on the integrated phys- 
ical spectral density or lattice data [Q]. We do believe it is the best way to confront the 
model under consideration with reality. 

The extraction of hadronic parameters from correlation functions by fitting is eas- 
ily beset with systematic errors, for example, because they are very sensitive to the 
parametrization of the fitting function. We have performed two types of fits: (i) us- 
ing the data points with statistical errors only (for results given in Table 1), and (ii) with 
statistical errors provided they are larger than 5 percent, otherwise they are put to be 5 
percent. The difference is related to a different weight assigned to the small x region, for 
which the statistical errors are very small. Generally, the difference in the fitted values is 
of the order of 10 percent, except in special cases like the p mesons, where the has a 
complicated valley in the parameter space, so that small changes in the calculated points 
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may drastically change the fitted parameters.. Lattice data typically resolve such problem 
by fitting only the long-distance part of the correlators, ignoring the small x behaviour. 

Having said this, we still would like to report the results of such fits. All points, from 
X = until the largest x, are fitted with 'pole plus continuum' expressions (two poles for 
axials), with only statistical errors included. A subset of results for diagonal 'fiavored' 
correlators is listed in Table 1. Let us make the following comments. 

- The masses of the lowest resonances in these channels, ranging from light pions to 
heavy axials, are very well reproduced by the model. The largest deviation, for the 
p meson, is about 15 percent. 

- The coupling constants A are somewhat smaller than the experimental values. It is 
probably due to the fact that the model does not have confinement, which effectively 
cuts off the tails of the wave functions and makes the particle more compact. The 
largest deviation is observed for the meson which is However, with our definition!™] 
of A, the model correctly reproduces the following natural trend: the coupling con- 
stants decrease from the 'most attractive' channels {a, vr) to the 'least attractive' 
ones (vectors and axials). In the most 'repulsive' channels (fiavored scalars) the 
resonances are absent altogether. 

- In contrast to the parameters of the lowest excitations, the position of 'continuum 
threshold' is surprisingly channel independent. Unfortunately, we do not have much 
phenomenological information about these important parameters (except for vector 
channels). However, if we assume that the threshold coincides with the position of 
the next to lowest state, the primed resonance, these numbers can be compared with 
the experimentally known 'primed' resonances. Indeed, in all known cases they are 
in the region around 1200 MeV, as the fitted thresholds indicate. Thus, one may 
speculate that the RIV model does predict the correct position of the 'continuum 
threshold' as well. 

^°Let us remind the reader, that it differs from standard definition in hterature. This allows us to 
compare different channels in a more straightforward way. 
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4.3. Summary 



We have reported a detailed calculation of a wide set of mesonic correlation functions 
in the 'random instanton vacuum'. The original idea was to use a crude model as a 
benchmark for further studies of interacting instantons, with the parameters (the instan- 



ton density and a typical radius) as suggested a decade ago ||10[, without any further 
adjustments. 

However, the results obtained show that this model works amazingly well, in partic- 
ular for the pseudoscalar correlators. Asymptotic freedom is recovered at small x, and 
at intermediated distances the correlation function can be parametrized by a spectral 
function with one pole (two for axials) and a continuum contribution. At the moment it 
is not clear why the phenomenological correlators are reproduced so well in spite of the 
apparent absence of such usual ingredients as confining forces and perturbative effects. 

In fact, not only the qualitative behaviour of all correlators at distances as large as 2 
fm is reproduced, but in a few cases, when a quantitative analysis is possible, we found 
agreement with the physical parameters at the 10-20 percent level. Thus, it is certainly a 
useful tool, so wider use of the model seems to be justified. 

At the same time, in two channels with a strong repulsive interaction, the rj' and 
the 6, the results are not satisfactory: the model correctly predicts such repulsion but 
overestimates its magnitude. We are looking forward to more realistic calculations with 
'interacting instantons' and hope to improve on this situation. 
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channel 


m^es [Mey] 


^/X[MeV] 


Eo [MeV] 


comment 


7r,0- 


142 ± 14 


510±20 


1360 ± 100 


this work 




138 


480 




phenomenology 




482 ± 12 


467±20 


1350±50 


this work 




495 




595 


phenomenology 


r/,0- 


500 
548 .8 


462 


1250 


this work 

phenomenology 




543 


500 


1160 


this work 




950±100 


390±20 


1500±100 


this work 




780 


409± 5 




phenomenology 


K\l- 


860±15 


341±20 


1300±50 


this work 




892 


448±25 




phenomenology 


0,1- 


850±50 


280±20 


1000±40 


this work 




1020 


492±15 




phenomenology 


ai, 1+ 


1132±50 


305±20 


1100±50 


this work 




1260 


400 




phenomenology 


K,,l+ 


1170±50 

1270 


302±20 


1120±50 


this work 

phenomenology 




1210±50 
1285 


293±20 


1200±50 


this work 
phenomenology 



Table 1 



The mass mres, the coupling constant A and the threshold energy Eq of some mesons in 
absolute units fitted according to the three-parameter formula for the spectral function. 
The pion mass is adjusted to the physical quark masses under the assumption that it 
is proportional to (m^ + md)^/^. Whenever available the parameters are compared with 
experimental numbers. 
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Figure Captions 



Fig. 1. Pseudoscalar correlation functions with quantum numbers of the pion and the 
kaon normalized to free massless quark correlator in the same channel plotted versus the 
distance x in fm. Our results are represented by the points. The phenomenological curves 
(as derived in [|l|) are given by the two solid lines. The dashed line corresponds to the 
'vacuum dominance' estimate (see text). Two dash-dotted lines show the fit described in 
the text. 

Fig. 2. Correlators for vector channels with quantum numbers of p, K* and (j) mesons 
normalized to the free massless quark correlator in the same channel plotted versus the 
distance x, in fm. Our results shown by squares, triangles and hexagons, respectively, 
should be compared with phenomenological expectations ||T[| (solid curves). The dashed 
line corresponds to the 'vacuum dominance' estimate. 

Fig. 3. Correlators for axial channels normalized to the free massless quark correlator in 
the same channel plotted versus the distance x in fm. Squares, triangles and hexagons 
show our results for the ud, us and ss channels, in this order. The first one should be 
compared with the phenomenological expectation given by the region in between the two 
solid curves The 'vacuum dominance' estimate is represented by the dashed curve. 
Three dash-dotted lines show the fits described in the text. 

Fig. 4. Correlators for isovector scalar channels normalized to the free massless quark 
correlator in the same channel plotted versus the distance x in fm. Squares and triangles 
show our results for the ud and us currents, respectively. The 'lower limit' to the phe- 
nomenological expectation (see text) is shown by the solid curve, while the dashed one 
corresponds to the 'vacuum dominance' estimate. The dash-dotted lines just guide the 
eye. 

Fig. 5. Correlators of the tensor current j^^ normalized to the free pseudoscalar massless 
quark correlator plotted versus the distance x in fm. Squares, triangles and hexagons 
show our results for the ud, us and ss channels, respectively. The 'vacuum dominance' 
estimate is represented by the dashed curve, and the dash-dotted curve shows the fit 
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described in the text. 



Fig. 6. The vector-tensor (a) and the pseudoscalar-axial (b) correlators in units of the 
pseudoscalar free quark correlator plotted versus the distance x in fm. The parameters 



of the solid curve in (a) are those obtained from QCD sum rules |14]. The dash-dotted 
line in (b) shows the fit described in the text. 

Fig. 7. Correlators for the isoscalar pseudoscalar channels t] and 77' normalized to free 
pseudoscalar massless quark propagator plotted versus distance x in fm. In figure (a) we 
show the correlators of the SU{3) singlet and octet currents (open and closed squares), 
while (b) represents the singlet-octet mixed correlator. Phenomenological expectations 
(see [1]) are shown by the solid curves. The dash-dotted curve in (a) is our fit for rj, and 
the dashed one shows the contribution of the r]' to the singlet correlator according to [1] . 

Fig. 8. Correlators for the isoscalar scalar a channel normalized to the free massless quark 
correlator in the same channel plotted versus the distance x in fm. The full correlator 
is shown in Figure (a), whereas (b) only has the connected part. The solid curve in 
(a) is the phenomenological expectation for the disconnected part. The long and short- 
dashed curves in (b) show the total correlator and the a— meson contribution to the fit, 
respectively. 

Fig. 9. Comparison between our results for the correlation functions (triangles) and the 
first lattice results [^] (squares). The four channels considered correspond to the isovector 
{e.g. ud) pseudoscalar, vector, scalar and axial the cases. The phenomenological result 
(solid lines) were derived in [|I|. All correlators are normalized with respect to the free 
quark correlator in the same channel and the distance on the a;— axis is in fm. 
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